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The prospect of a Dirac half metal, a material which is characterized by a bandstructure with a 
gap in one spin channel but a Dirac cone in the other, is of both fundamental interest and a natural 
candidate for use in spin-polarized current applications. However, while the possibility of such a 
material has been reported based on model calculations[H. Ishizuka and Y. Motome, Phys. Rev. 
Lett. 109 , 237207 (2012)], it remains unclear what material system might realize such an exotic 
state. Using first-principles calculations, we show that the experimentally accessible Mn intercalated 
epitaxial graphene on SiC(OOOl) transits to a Dirac half metal when the coverage is > 1/3 monolayer. 
This transition results from an orbital-selective breaking of quasi-2D inversion symmetry, leading to 
symmetry breaking in a single spin channel which is robust against randomness in the distribution of 
Mn intercalates. Furthermore, the inclusion of spin-orbit interaction naturally drives the system into 
the quantum anomalous Hall (QAH) state. Our results thus not only demonstrate the practicality of 
realizing the Dirac half metal beyond a toy model but also open up a new avenue to the realization 
of the QAH effect. 

PACS numbers: 75.50.Pp, 73.43.-f, 81.05.ue 

The discoveries of novel band structure^ such as 
the Dirac spectrum [ll-[^ and half-metallicity[l, Q, have 
sparked tremendous scientific and technological inter¬ 
est. By combining the two fascinating properties, a po¬ 
tentially more interesting state, namely, the Dirac half 
metal, has been proposed recently based on a model 
calculation. 0 Such a system is not only semiconduct¬ 
ing in one spin channel and metallic in the other but 
also characterized by the zero-gap band structure with 
a linear dispersion. The coexistence of 100% spin polar¬ 
ization and massless Dirac fermions makes this system 
promising for future applications in electronics, spintron- 
ics and optoelectronics. In addition, as the Dirac half 
metal possesses a Dirac cone in only one spin channel 
(being an ordinary insulator in the other spin-channel), 
the gap opening triggered by spin-orbit coupling (SOC) 
leads to topological phase transition in only one spin 
channel. Hence, the Dirac half metal is also a natural 
avenue toward the realization of the quantum anomalous 
Hall effect (QAHE). 

Symmetry plays a central role in Dirac materials, such 
as the sublattice symmetry for graphene and the time- 
reversal symmetry for topological insulators. [iSIl It pro¬ 
tects two linear dispersive bands crossing at the Dirac 
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point and its breaking leads to a Dirac gap, which is 
an essential requirement for modern electronics. From 
this symmetry point of view, the Dirac half metal may 
be more interesting because there must exist a symmetry 
order that is destroyed in only one spin channel but not in 
the other, which thus calls for a distinctly different sym¬ 
metry property from that of the non-spin-polarized Dirac 
materials like graphene and topological insulators. Given 
the great academic interest and potential applications, it 
is highly desirable to search for experimentally realiz¬ 
able Dirac half metals with simple crystal and electronic 
structures. As the half metal has broken time-reversal 
symmetry alongside with a spin-resolved orbital physics, 
a viable strategy for designing the Dirac half metals is 
to search for such an orbital physics, whereby the origi¬ 
nally degenerate orbitals, e.g., the transition metal (TM) 
d-orbitals, split into two groups with opposite behaviors 
under the same symmetry operation. 

In this work, we introduce the first Dirac half metal 
system to be engineered through such a mechanism, Mn- 
intercalated graphene on SiC(OOOl). By utilizing the sub¬ 
strate modulation, the characteristic Mn d-orbitals hy¬ 
bridizing with graphene p-orbitals in the two spin chan¬ 
nels become different, so that the corresponding p-d inter¬ 
action energies acquire, respectively, odd and even par¬ 
ities under the same quasi-2D inversion operation. Fur¬ 
thermore, as the substrate orders a subset of the Mn 
atoms, the symmetry-breaking induced gap-opening of 
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the majority spin channel is global and the symmetry 
cannot be recovered through a random distribution or 
local fluctuations of Mn atoms. Using first-principles 
calculations, we show that this system becomes a Dirac 
half metal when the Mn coverage y exceeds 1/3 mono- 
layer (ML: Note that 1 ML coverage is defined as one 
Mn atom per surface Si atom.), having zero gap in the 
minority spin channel while exhibiting a 150-meV gap in 
the majority spin channel. The experimental realization 
of such a system is supported not only by our calcula¬ 
tions which indicate that the Dirac half metal is robust 
against Mn clustering and geometric distortion, but also 
by previous work which has demonstrated the experimen¬ 
tal feasability of Mn intercalation into epitaxial graphene 
on SiC(0001)[l3. Finally, the inclusion of SOC in the 
calculation reveals that the Dirac half metal exhibits the 
QAHE, as supported by our Chern number analysis. 


Spin-polarized calculations were performed using the 
Vienna ab initio simulation package (VASP) [l^ within 
the framework of the density-functional theory (DFT). 
The local density approximation (LDA)[I| and the 
projector-augmented wave [l^ potential with a cutoff 
energy of 400 eV were used to describe the exchange- 
correlation energy and the electron-ion interaction, re¬ 
spectively. A fc-mesh of 6 x 6 x 1 was used to sample the 
Brillouin zone. The adopted model was to intercalate dif¬ 
ferent coverages of Mn atoms into the epitaxial graphene 
on SiC(OOOl), which was based on a x •\/3i?30° unit 
cell for 6iL-SiC to accommodate a 2 x 2 graphene. In the 
simulations, the bottommost three out of the six SiC bi¬ 
layers were fixed at their respective bulk positions while 
all other atoms were fully relaxed without any symmetry 
constraint until the residual forces were less than 0.01 

eV/A. 

Quasi-2D inversion symmetry.-The quasi-2D inversion 
operation is given as follows: —— and 

where O is the inversion center of the graphene, 
lies in the plane of graphene lattice sites and 
is normal to the plane (see Fig. 2). A physical quantity 
possesses the quasi-2D inversion symmetry, only if it is 
invariant under such an operation. As the graphene Dirac 
cone is related to its inversion symmetry, introducing the 
quasi-2D symmetry helps to understand the effect of the 
TM on the graphene Dirac cone. The Mn-intercalated 
graphene on SiC(OOOl) is characterized by the p-d hy¬ 
bridization between the Pz oribtals of graphene and the 
d oribtals of the TM = Mn. Within the Slater-Koster 
approximation [l6j . the p-d interactions exhibit two kinds 
of symmetry properties, because, upon the quasi-2D in¬ 
version, only the in-plane, x and y-related directional 
cosines change signs, while the out-of-plane ^-related di¬ 
rectional cosines do not. This leads to a sign change in 
the Pz (Cydxz (TM) andpz (C)-dy 2 (TM) interaction en¬ 
ergies, while the signs for the Pz (C)-dxy, Pz iC)-dx 2 -y 2 , 
and Pz {C)-dz 2 interaction energies will not change. In 



FIG. 1: (Color online)Illustration of the quasi-2D inversion 
operation to characterize the symmetry physics for a transi¬ 
tion metal (Green ball) atom adsorption above the graphene 
hexagonal center. Under the quasi-2D inversion operation 
{R —^ —R), the sign changes only for x and y vectors while 
not for z vector. So do their directional cosines ri {i = x, y 
and z). Within the Slater-Koster approximation, the p-d in¬ 
teraction energies are (pz\H\dxy) = rxVyrzi'/iVpdc - 
{Pz\H\d^2_y2) = \rz{rl-rl){y/Wpda-2Vpdz), {pz\H\d^2) = 
rz[{rl - (rl + rl)l2)Vpda + y/iirl + ry)V'pd^], (j)z\H\dxz) = 
rxWirlVpda + {t-2rl)Vpd-K\, and {pz\H\dyz) = ry[s/irlVpda + 
(1 — 2r1)VpdTT\, respectively, where Vpda- and Vpdz are two pa¬ 
rameters. These clearly show that the sign change occurs un¬ 
der the quasi-2D inversion operation for pz hybridization with 
either dxz or dyz orbitals, meaning the symmetry breaking. 
In contrast, they are unchanged for the other three d-orbitals 
{dz'z, dxy and d^ 2 _y 2 ), meaning the symmetry invariant p-d 
interaction. 


other words, the dxz (TM) and dj, 2 (TM) orbital inter¬ 
action with graphene Pz orbitals leads to a perturbation 
which breaks the inversion symmetry graphene, leading 
to a Dirac gap.This is not true for the other d-orbitals. 
In this regard, a Dirac half metal could be realized in 
our system if the d-orbitals in each spin channel that 
participate in the p-d hybridization can be adequately 
engineered. 

Dirac half metal behavior in the Mn intercalated 
graphene.- 

The X = 1/3 ML Mn coverage, corresponding to a 
single Mn for every three Si at the SiC(OOOl) surface, 
is a critical coverage in this system. At x = 1/3 ML, 
each Mn is anchored by three Si and all of the Si dan¬ 
gling bon ds ( DBs) are passivated and the system is non¬ 
magnetic. [l7|,ll^ Further insertion of Mn inevitably leads 
to the formation of Mn trimmers, as illustrated in Figs. 
2(a) and 2(b), where the optimized structure for x = 
5/12 ML is shown. The 5/12 ML coverage is equivalent 
to adding a fifth Mn atom (Mns) to a 2 x 2 supercell 
with X = 1/3 ML. 

Figures 2(c) and 2(d) show the spin-resolved band 
structures for x = 5/12 ML. A notable feature is the fully 
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FIG. 2: (Color online)Geometric and electronic structures of 
Mn intercalated epitaxial graphene on SiC(OOOl). (a) Top and 

(b) side views of the optimized geometry for Mn intercalation 
coverage of y; = 5/12 ML. Blue rhombus in (a) denotes the 
supercell. The numbers (1, 2, 3, 4, and 5) mark the five 
different Mn atoms. Note that there is only one configuration 
for adding Mns into a 2 x 2 supercell with x = 1/3 ML. Here 
Mns, Mn 4 and Mns atoms form a trimer. The corresponding 
spin-resolved band structure along the high-symmetry lines: 

(c) majority and (d) minority spin channels, (e) Total density 
of states as well as the partial contributions from graphene 
(G), intercalated Mn and surface Si atoms. The Fermi level 
is set to zero (brown dot). Note that spin-orbit coupling is 
neglected in this calculation. 


spin-polarized Dirac spectrum. The majority spin chan¬ 
nel has an energy gap of 190 meV, while the minority spin 
channel has a Dirac cone off from K. A displacement of 
the Dirac point to a lower-symmetry point has been fre¬ 
quently seen in the graphene-related systems, e.g., when 
an anisotropic strain is present. IT], [l^ To confirm this 
is indeed a Dirac point, we have calculated the three- 
dimensional band structure (not presented here) to find 
that the level crossing does happen at the Fermi level. 
Moreover, by using a 2 x 1 supercell, we find the ferro¬ 
magnetic coupling configuration is about 15 meV more 
stable than the antiferromagnetic one, showing the long 
range ferromagnetic order. 

The Dirac half metal is also revealed by the density of 
states (DOS) in Fig. 2(e). Here, the majority spin gap is 
150 meV, which is smaller than the 190 meV at K in Fig. 


2(c), suggesting that the minimum gap may also be off 
K. The projected DOS reveals that the states around the 
Fermi level are dominated by Mn, indicating its essential 
role in forming the half-metallicity and linear dispersion. 
Unlike conventional Dirac systems, here an excited Dirac 
fermion can be fully spin-polarized with unique advan¬ 
tages in polarization optics and spintronics.[20| 

To explain the first-principles results, let us consider 
the following two questions: (1) how does the system 
magnetism arise, and (2) why does the gap open only in 
the majority-spin channel? Answer to the first question 
lies in the different interactions of the Mn atoms with the 
SiC substrate, whereas the answer to the second question 
lies in the quasi-2D inversion symmetry discussed previ¬ 
ously. 

For the magnetism, we note that there are only 12 
surface Si atoms per unit cell at % = 5/12 ML, each 
with 1 DB. The stable form of Mn on SiC(OOOl) is to 
bind with 3 Si DBs in a r 4 conhguration, as we have 
demonstrated before 17|. This means at 1/3 ML, all the 
DBs have been saturated. The lack of sufficient Si DBs 
at 5/12 ML naturally divides the Mn atoms into two 
categories: (i) Mni-Mn 4 , each is chemically bound to 
three Si atoms and (ii) Mns, which is atom-like because 
of no more surface Si to bind to, in spite of that the 
crystal field effect promotes all its 7 valence electrons 
into the 3d orbitals. For category-(i) Mn, binding to 
Si completely quenches its spin-polarization[l3. Among 
the crystal-field-split d-orbitals, the z-involving orbitals, 
dxz, dyz, and dz 2 ^ are strongly hybridized with the Si 
DBs to result in significantly lower energies, as shown 
in Fig. 3(a). The remaining d^y and dx 2 _y 2 orbitals 
are, however, little affected by the Si from the substrate 
and hence they maintain their non-bonding characters 
with energies near the Dirac point of the graphene. For 
category-(ii) Mn, on the other hand, as none of its d- 
orbitals can be significantly affected by the substrate, 
they are all atom-like, i.e., non-bonding with energies 
near the Dirac point of the graphene, as shown in Fig. 
3(b). In this case, the 7 d-electrons should occupy 5 spin- 
up states and 2 spin-down states to result in a high spin 
configuration of 3 /is, as confirmed by our first-principles 
calculation. 

For the gap opening, we consider the coupling between 
Mn d and graphene Pz orbitals. For substantial coupling 
to exists, the states need to be close in energy, namely, 
between the non-bonding d orbitals highlighted in the 
green areas in Figs. 3(a) and 3(b) and pz orbitals in 
Fig. 3(c). For the spin-minority channel, the coupled 
states only involve dxy and d,j, 2 _y 2 , whose hybridization 
with Pz orbitals holds the quasi-2D inversion symmetry, 
which protects the Dirac cone. In contrast in the spin- 
majority channel, the symmetry-breaking interactions in¬ 
volving the dxz and dyz orbitals open a Dirac gap. Put 
together, this analysis shows that the Dirac half metal 
arises due to the occupation-modulation effect by SiC 

















4 


Si-anchored Mn(l,2,3,4) magnetic impurity Mns Graphene 


majority 

minority 

(b) 

majority minority 

d}cz. dy?-, dz- 

(c) 

V Dirac cone 

dxy, 

dy.y. 


daz, dvz. dz- 




Hr H" 


LLLLL + + 





d\y, dx--y- dxy. dx-y- 



e 

+ + + drdrdr 

dxz, dy?., dz- duz, dyz. dz- 




M 



FIG. 3: (Color online)The d-electron occupations of two cat¬ 
egories of Mn atoms and the schematic diagram of their elec¬ 
tronic levels with respect to graphene Dirac point, (a) Mni 
to Mn 4 , fully saturated by three Si DBs. (b) Mns, atomic-like 
magnetic impurity, (c) Graphene Dirac cone. The relative en¬ 
ergy alignments are deduced from the symmetry analysis in 
combination with our DFT calculations. Green box denotes 
those states important for the p-d coupling, (d) Schematic il¬ 
lustration of the disorder due to possibly random distribution 
of Mns. For clarity, only graphene and Mn (pink balls) are 
shown. Blue rhomboid denotes the 2x1 supercell containing 
two Mn trimers used for test calculations on random disorder, 
(e) Geometry at x = 1/3 ML, where there exist category (i) 
Mn atoms underneath the yellow hexagons. Because of this 
pre-established distribution of Mn, the symmetry breaking 
cannot be canceled on average and only the one on the sites 
(in triangle) bonding-free to category (i) Mn is anticipated to 
be suppressed by the random effect. 


substrate on the Mn d orbitals and the quasi-2D inversion 
symmetry of the combined graphene-plus-Mn system, as 
shown in Fig. 1. 

We calculate the Mn intercalation energy, which is de¬ 
fined as the total energy difference between (Mn bulk -|- 
epitaxial graphene on SiC) and the Mn intercalated sys¬ 
tem. We obtain 0.4, 0.3, and 0.0 eV/Mn, for y = 1/3, 
5/12, and 1/2 ML, respectively. The higher stability at 
low X suggests that Mn clustering is energetically unfa¬ 
vorable. In other words, Mn tends to avoid each other, 
which allows them to maximize the binding to available 
Si sites on the substrate. Thus, trimer formation primar¬ 
ily happens after all Si sites have been fully passivated, 
i.e., 1/3 ML. 

This means only when y > 1/3 ML, there emerges the 
category-(ii) Mn and the system becomes a Dirac half 
metal. As y gradually increases to 5/12 ML, more and 
more Mn trimers form, and the system remains a Dirac 
half metal according to our analysis based on the two 
categories of Mn. This is explicitly corroborated by our 
first-principles calculations on the y = 17/48 ML system 
(containing only one Mn trimer in a 4-\/3 x 4^37730° SiC 


lattice). Increasing the Mn coverage so that all of the Mn 
atoms form trimers (y = 1/2 ML), however, results in a 
system which is no longer a Dirac half metal, and instead 
exhibits charge transfer from the majority to minority 
spin channel and a non-integer total moment of 2.8 pLs- 
These calculations indicate that the Mn intercalated SiC 
(0001) system possesses a wide range of coverages (1/3 
ML< y < ymax), with (5/12 ML < ymax < 1/2 ML), in 
which it exists in the Dirac half metal state. 

The use of supercells in these calculations raises the im¬ 
portant question as to whether the effects of the symme¬ 
try breaking will diminish, or be completely eliminated, 
in realistic systems in which the Mn trimers may dis¬ 
tribute randomly, as schematically shown in Fig. 3(d). 
This is indeed the case for systems in which the symme¬ 
try is only locally broken, and the averaging of the per¬ 
turbations leads to a situation which preserves the bulk 
symmetry of the graphene on substrate 21 - 2^ . While 


our first-principles calculations suggest that the gap per¬ 
sists in the presence of disorder (e.g. the 2x1 supercell 
[blue rhomboid in Fig. 3(d)] with reduced order is not 
only ^13 meV less stable than the one in Fig. 2(a) but 
still yields a similar half metallic gap), this is not a ques¬ 
tion which can be answered, strictly speaking, based on a 
supercell approximation which imposes long range order. 
The important realization is that the substrate fixes the 
positions of category-(i) Mn [Mni-Mn 4 ] and as a result 
the symmetry breaking perturbation of category-(ii) Mn 
[Mns] cannot be completely random. 

To illustrate this last point, consider the worst case sce¬ 
nario for maintaining a symmetry breaking perturbation, 
namely, complete randomness of the distribution of Mns 
atoms on all of the possible Mns sites (centered on the 
white hexagons in Fig. 3(e)). If one adds all of the asym¬ 
metric Slater-Koster interaction terms effecting graphene 
resulting from such a distribution, one finds that they 
completely cancel only on the C-atoms marked with a tri¬ 
angle in Fig. 3(e). Note that all of the C-atoms adjacent 
to a category-(i) Mn [Mni-Mn 4 , centered on the yellow 
hexagons] are still effected by symmetry breaking per¬ 
turbations. The reason complete cancellation does not 
occur is because Mns cannot occupy the yellow hexagon 
sites. This regular lattice of sites is determined by the 
substrate and any Mn which sits at one of these sites is of 
category-(i) due to its interaction with the substrate. As 
a result, even a completely “random” distribution of Mns 
still has long range order in that there is a regular lattice 
of sites from which it is excluded. Hence, the symmetry 
breaking of this system is inherently global and cannot 
be eliminated due to randomness in the distribution of 
Mns. 

Quantum anomalous Hall effect.-lt has been shown 
previously that Mn intercalation significantly enhances 
the SOC. At X = 1/3 ML, it opens a band gap of 26 
meV. This drives the system into a quantum spin Hall 
state. 18| In the present case, Mns makes the system spin- 
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FIG. 4: (Color online) Band structure including the SOC 
(top panel) and the Berry curvature (bottom panel) along 
the high-symmetry line of Mn intercalated epitajdal graphene 
on SiC(OOOl) at X = 5/12 ML. Red circles highlight the SOC 
gap. The Fermi level is set to zero. 


polarized, thereby breaking the time reversal symmetry. 
In this case, the SOC should open the Dirac cone in the 
minority spin channel with a nontrivial gap and thus the 
system enters a QAH state. First-principles calculations 
including the SOC indeed yields a SOC gap around 20 
meV at X = 5/12 ML, as shown in the top panel of 
Fig. 4 (Red circles). To identify the topological prop¬ 
erties of the gapped state, we employ the Chern num¬ 
ber analysis [2^ [26| in which the Chern number is calcu¬ 
lated by a fc-space integral over the first Brillouin zone, 
C = -^ /bz ^(k)d^k, where r2(k) is the Berry curvature 
of all occupied states [27|: 


= - E E 

ndEp m^n 


(V^nk l^a: I'^rnk) ('^rrik | | V^nk) 

(S'mk '^nk)^ 


( 1 ) 

ipnk is the spinor Bloch wave function of band n with cor¬ 
responding eigenenergy e„k, and Vi is the zth Cartesian 
component of the velocity operator. The calculations 
are conducted by using the WANNIER90 package 2^ 281 
with the maximally localized Wannier functions, which 
reproduces the DFT results well (not shown). Figure 
4 (bottom panel) shows the calculated total Berry cur¬ 
vature along the high symmetry lines. The large peaks 
between F and K {—K) points arise where the intersec¬ 
tions between the conduction and valance bands split by 
the SOC, giving rise to small denominators in Eq. (1), 
and hence a large contribution to r2(k). By integration 
over the Brillouin zone, we obtain an odd Chern number, 
C = 1, to confirm that, in a ribbon structure, only one 
chiral state can exist on the edge. It should be noted 
that this Chern number is different from that in TM- 


doped graphene [ 29 ! and silicene[3lj| where the Chern 
number is even ((7 = 2). 

To summarize, we propose a strategy to make Dirac 
half metals in experimentally-accessible Mn intercalated 
epitaxial graphene on SiC(OOOl), whose low-energy ex¬ 
citations have a spin-gapless Dirac spectrum, originated 
from the coupling between the substrate-modulated Mn 
d orbitals and graphene Pz orbitals. In the majority spin 
channel, the interactions involving Mn dxz and dyz or¬ 
bitals break the quasi-2D inversion symmetry and thus 
induce a Dirac gap. In the minority spin channel, how¬ 
ever, there exist only dxy and dx 2 -y 2 interactions, which 
leaves the Dirac cone intact. Furthermore, the substan¬ 
tial SOC of Mn drives the Dirac half metal system into 
a QAH state. As the SOC is an intrinsic property of 
the material, the QAHE is also an intrinsic property of 
the Dirac half metal. We thus anticipate that our study 
opens a new route toward the realization of the QAHE. 
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